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ABSTRACT 

We present multiwavelength observations of the persistent Fermi-LAT unidentified 7 -ray source 
IFGL J1417.7—4407, showing it is likely to be associated with a newly discovered X-ray binary 
containing a massive neutron star (nearly 2Mq) and a ^ 0.35Mq giant secondary with a 5.4 day 
period. SOAR optical spectroscopy at a range of orbital phases reveals variable double-peaked Ha 
emission, consistent with the presence of an accretion disk. The lack of radio emission and evidence 
for a disk suggests the 7 -ray emission is unlikely to originate in a pulsar magnetosphere, but could 
instead be associated with a pulsar wind, relativistic jet, or could be due to synchrotron self-Compton 
at the disk-magnetosphere boundary. Assuming a wind or jet, the high ratio of 7 -ray to X-ray 
luminosity 20) suggests efficient production of 7 -rays, perhaps due to the giant companion. The 
system appears to be a low-mass X-ray binary that has not yet completed the pulsar recycling process. 

This system is a good candidate to monitor for a future transition between accretion-powered and 
rotational-powered states, but in the context of a giant secondary. 

Subject headings: pulsars: general — Gamma rays: general — X-rays: general — binaries: spectro¬ 
scopic 


1. INTRODUCTION 

Millisecond pulsars are thought to have been recycled 
to short periods through accretion from either a main se¬ 
quence or giant secondary. Despite the expectation that 
this mass transfer process should be lengthy, most known 
field millisecond pulsars have degenerate white dwarf 
companions that represent the end stage of the recycling 
process (see the review of Tauris & van den Heuvel 2006). 
This suggests the radio emission is quenched during even 
low-level accretion onto the neutron star, making the sys¬ 
tems undetectable as pulsars. 

The discovery of PSR J1023-I-0038 represented the first 
compelling evidence that some millisecond pulsars ac¬ 
tively switch between a low state with radio pulsations 
and a high state with an accretion disk (Archibald et 
al. 2009), which in this system is due to Roche lobe over¬ 
flow of its ~ 0.2Mq bloated main sequence companion 
(Archibald et al. 2013; Bogdanov et al. 2015). Two ad¬ 
ditional similar systems have been discovered, both with 
main sequence companions (Papitto et al. 2013; Bassa et 
al. 2014; Roy et al. 2015). 

These transitional systems belong to a larger class of 
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objects known as “redbacks”, binary systems in which 
the wind of the pulsar primary ablates a main sequence 
companion of > 0.1 Mq, causing radio eclipses for a sub¬ 
stantial fraction of the orbit. The Fermi Large Area 
Telescope (LAT) has been a boon for the discovery of 
redbacks: many of these are 7 -ray sources that had not 
previously been detected as pulsars, possibly because of 
the presence of ionized material in the system (Roberts 
2013). 

As part of an ongoing survey of unidentified Fermi- 
LAT sources (Strader et al. 2014), this paper reports ob¬ 
servations of a faint 7 -ray source from the IFGL catalog: 
IFGL J1417.7-4407 (Abdo et al. 2010). Using X-ray, ra¬ 
dio, and optical observations including photometry and 
spectroscopy, we find that IFGL J1417.7—4407 is likely 
a Galactic compact binary with a massive neutron star 
primary and a giant secondary in a relatively wide orbit. 

2. OBSERVATIONS 
2 . 1 . Fermi-LAT Source 

The 7 -ray source was first catalogued as IFGL 
J1417.7—4407 (|Abdo et al.ll2010D . detected in 11 months 

of LAT all-sky survey dat a. _ In the recently released 

3FGL catalog (|Acero et al.ll2015h . deriyed from 4 years 
of LAT obseryations, it is listed as 3FGL J1417.5—4402. 
The 3FGL 95% error ellipse (3.7' x 3.5' at position 
angle 17°), centered at epoch J2000 (R.A., Dec.) = 
(214.377°, —44.043°), is 9x smaller in area than 
the IFGL error ellipse. The 3FGL source is character¬ 
ized by a single power-law spectrum with an index of 
2.37 ± 0.08, with no significant eyidence for a curyed 
spectrum. The 0.1-100 GeV flux corresponds to a lu¬ 
minosity of (2.8 ± 0.3) X 10^^ (fi/4.4kpc)^ erg s“*. The 
3FGL catalog light curye appears superhcially yariable, 
but the yariability index is 55 (see Nolan et al. 2012 for a 
definition), below the 99% threshold of ^ 72 necessary to 






2 


Strader et al. 


classify the variability as “probable” in the 3FGL cata¬ 
log. Therefore we assess the 7 -ray variability as uncertain 
with present data (though there is evidence for long-term 
X-ray variability; see below). 

2 . 2 . X-ray Observations 

IFGL J1417.7—4407 was observed with Chandra as 
part of a Gycle 12 program (PI; C. Ricci) targeting six 
unidentified IFGL Fermi-LAT sources. The ~ 2 ksec ex¬ 
posure was obtained on 2011 Feb 10, and had placed the 
target on the ACIS-I detector (~ 17' x 17' field of view) 
to ensure the entire IFGL J1417.7—4407 error ellipse 
was covered. The XASSIST processed catalog source, 
X141730.55—440257.6 (J2000 position based name; see 
iPtak fc Griffiths! I2003H . was the only X-ray source de¬ 
tected in the Chandra exposure within the 3FGL Fermi- 
LAT 7 -ray error ellipse (§ 2.1), being located essentially 
at its centroid (^ 0.4' offset). 

We downloaded the Chandra data from the archive 
and analyzed the data using the GIAO software v.4.7 
(|Fruscione et al.l 12006(1 and the most recent CALDB 
v4.6.5. The level 1 (evtl) files were reprocessed using 
standard procedures in CIAO to generate new evt2 files. 
We found no background flares during the observation, 
thus the net exposure was 2051.2 s from 09:43:10.032 
to 10:17:21.232 UTG (</> = 0.94 using the ephemeris from 
§3). For spectral analysis, we extracted photons with en¬ 
ergies 0.5 — 7 keV within a circle with radius 5" centered 
on the X-ray position, with the background determined 
from an annulus from 15" to 160", giving a net count 
rate of 26.1 ± 3.6 cts ksec“^ (53 photons). 

We found that an absorbed single power-law fit to the 
data could not be distinguished from an absorbed black- 
body spectrum (cstat/dof = 192.7/444 and 192.2/444, 
respectively). The best-fit power-law had a photon in¬ 
dex F = 1.32 ± 0.40 and an observed 0.5 — 7 keV flux 
of (5.6 ± 1.4) X 10“^^ erg cm“^ s“^ (90% confidence), 
while the blackbody had kT = 0.781°^eV with flux 
(4.8llg) ^ 10“^° erg cm“^ s“^. These values were de¬ 
rived a ssuming the Galactic absorption, Nh = 6.54 x 10^° 
cm“^ (|Kalberla et al.ll2005ll . and gave unabsorbed fluxes 
5.9 X 10“^° and 5.0 x 10“^° erg cm“^ s“^ for the respec¬ 
tive models. There were hints of additional absorption 
but when left to vary, the errors are compatible with 
zero, thus only the Galactic absorption fixed values are 
presented. The power-law flux is equivalent to a lumi¬ 
nosity of (1.4±0.4) X 10°°((i/4.4kpc)^ erg s“^, where the 
uncertainty is at 90% confidence as above. 

To gauge possible variability, we divided the counts 
into two bins of 1025.6 s each and found net count rates of 
36.5±6.0 (37 photons) and 15.7±4.0 (16 photons) counts 
ksec“^ in the respective bins. Assuming the power-law 
model, the unabsorbed 0.5 —7 keV fluxes are (6.9±1.9) x 
10 “^° and 4 .OII 7 ^ 10“^° erg cm“^ s“^, respectively, and 
are consistent within the joint uncertainties. Due to the 
low count rate and unfortunate location of the source 
over a node in an AGIS-I chip, these data are not well- 
suited to provide constraints on the shorter-term X-ray 
variability observed for known transitional millisecond 
pulsars (e.g., Bogdanov et al. 2015). 

We note that ROSAT catalogued an X-ray source, 
2RXP J141731.0-440253 (observation on 1997 Feb 09; 
IROSAT Consortium! I 2 OOOII that is only 6.5" offset from 


the Chandra source position, thus well withi n the typical 
posit ional uncertainty for ROSAT sources (iVoges et al.l 
II 999 II . The ROSAT PSPC count rate of this source is 
42.4±8.0 counts ksec“^ over 0.1-2.4 keV, which is equiv¬ 
alent to observed 0.5-7 keV fluxes of 1.7 x 10“^^ erg 
cm“^ s“^ and 2.0 x 10“^^ erg cm“^ s“^ assuming the 
best power-law and blackbody models derived from the 
Chandra data. These are a factor of 3 or more brighter 
than the corresponding Chandra fluxes, thus there is ev¬ 
idence for long-term X-ray variability of the source. The 
source has not been detected by all-sky X-ray monitors 
such as Swift/BAT (Baumgartner et al. 2013) or MAXI 
(Hiroi et al. 2013), suggesting it has not undergone an 
outburst in the recent era of these missions. 

2.3. Radio Observations 

Petrov et al. (2013) presented radio continuum obser¬ 
vations at 5.5 and 9.0 GHz with the Australian Telescope 
Compact Array of a number of unidentified Fermi-LAT 
sources. During these observations, obtained on 2012 Sep 
20, no radio emission was detected at the position of the 
X-ray and optical counterpart discussed in this paper. 
The only radio source found within the 3FGL ellipse is 
a faint flat-spectrum radio source (2.6 ± 0.3 mJy at 5.5 
GHz; a = -kO.2; oc :/“) at R.A. = 14:17:16.39, Dec. = 
—44:04:39.4 offset by ^ 3.1' from the X-ray/optical po¬ 
sition. 

This field was re-observed using the 6 km configura¬ 
tion on 2013 Feb 4, with two 15.5 min pointings simul¬ 
taneously observed again at 5.5 and 9 GHz, each with a 
bandwidth of 2 GHz. Again, no radio emission associated 
with the X-ray/optical source was detected, with a flux 
density upper limit of 0.15 mJy at both frequencies. The 
field radio source was found to be significantly variable 
with a 5.5 GHz flux density of 1.43±0.07 mJy (a = -0.2) 
in the 2013 observation; it is not detected in the Chandra 
data with a 0.5-7 keV X-ray upper limit of < 1.5 x 10“^"^ 
erg cm“^ s“^. The results of Schinzel et al. (2015) imply 
a probability of 15% to find an unassociated 2 mJy com¬ 
pact radio source within 3.1' of a Fermz-LAT source. In 
addition, extragalactic high-latitude Fermi-LAT sources 
are mostly blazars with radio flux densities 10-100 times 
higher than the 2 mJy measured for this source (Acker- 
mann et al. 2015). Thus it is reasonable to conclude this 
radio source is unrelated to the 7 -ray source. 

2.4. The Optical and Near-IR Counterpart 

Using the USNO Bl.O catalog (jMonet et al.ll2003l l. we 
found an optical counterpart within 0.63" of the X-ray 
source, located at a J2000 sexigesimal position of (R.A., 
Dec.) = (14:17:30.604, -44:02:57.37). This source has 
optical magnitudes B2 = 16.90, R2 = 15.79, and / = 
14.98 mag. It is coincident with a 2MASS point source 
with J = 14.17 ± 0.03, H = 13.50 ± 0.03, and K = 
13.41 ± 0.04 (Cutri et al. 2003). It is also present in the 
UCAC4 astrometric catalog (Zacharias et al. 2013), with 
a {fiaCOsS, fis) proper motion of (- 8.8 ± 5.0, -3.1 ± 4.7) 
mas yr“^. 

2.4.1. Catalina Sky Survey 

We obtained archived photometric observations of the 
USNO Bl.O sourc e from the Gatalina Sky Survey (CSS; 
(Drake et al.l (200911 and its associated southern Siding 
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Spring Survey (SSS). These surveys take unfiltered im¬ 
ages in sequences of four 30-sec exposures typically reach¬ 
ing V"-equivalent magnitudes of ~ 19-20 mag. We found 
the source catalogued in SSS as J141730.6-440257, with 
195 photometric measurements from 2005 Aug 15 to 2013 
May 7. Four of these had large uncertainties (> 0.15 
mag) and were removed, leaving a total of 191 measure¬ 
ments. The median magnitude is t4quiv = 15.85 and 
median uncertainty 0.09 mag. There is no evidence for 
a change in the mean brightness of the source over the 
range covered by these data. 

Using the Lomb-Scargle periodogram as implemented 
in R (Ruf 1999), we searched for the most likely pho¬ 
tometric period in the SSS photometry. This period is 
2.69 d, but it turns out that this is an alias of the real 
period determined via spectroscopy (§3.1). This period 
is twice as long: 5.3737 ± 0.0003 d, with the uncertainty 
estimated via bootstrap. 

2 . 4 . 2 . PROMPT 

We obtained time series photometry of the optical 
source in BVR with the PROMPT-5 telescope (Reichart 
et al. 2005) at Cerro Tololo International Observatory be¬ 
tween 2014 June 15 and 2014 August 29. The data were 
reduced in the standard manner, and aperture photome¬ 
try performed to obtain instrumental magnitudes for the 
target source differentially with respect to five compari¬ 
son stars in the field. These magnitudes were calibrated 
using observations of Landolt (1992) standard star fields 
observed over a range of airmasses on most nights, with 
zeropoints and extinction coefficients determined on pho¬ 
tometric nights. 

2 . 5 . Optical Spectroscopy 
2 . 5 . 1 . SOAR Spectroscopic Monitoring 

We obtained two low-resolution spectra of the candi¬ 
date optical counterpart to IFGL J1417.7—4407 on 2013 
Jan 15 using the Goodman High-Throughput Spectro¬ 
graph (Glemens et al. 2004) on the SOAR 4.1-m telescope 
using a 1.03" slit and a 600 1 mm“^ grating (resolution 3 
A) . The spectrum appears to be that of a late G or early 
K star (see Figure 1), with the exception of a bright Ho? 
emission line with an equivalent width of ~ 8 A. The 
Ha emission line has a double-peaked morphology, with 
the centroids of the blue and red peaks each offset by 
^ 120 km s“^ with respect to the photospheric velocity 
of the star. This emission is consistent with the presence 
of an accretion disk, and considering its X-ray detection, 
we concluded this source was likely to be the true optical 
counterpart of the Fermi source. 

The following season we began spectroscopic moni¬ 
toring of this optical counterpart, taking spectra in 14 
epochs from 2014 Jan 11 to 2015 Feb 18. Most obser¬ 
vations were taken with a 1.03" slit and a 1200 1 mm“^ 
grating centered at 6100 A, yielding a resolution of 1.7 A 
and wavelength coverage of ~ 5375-6640 A. This enabled 
measurements of the photospheric radial velocity as well 
as monitoring of the Ha emission. At a few epochs, due 
to the unavailability of the 1200 1 mm“^ grating, we in¬ 
stead used a 2100 1 mm“^ grating, which had a reso¬ 
lution of 0.9 A and wavelength coverage ~ 4960-5600 
A. The exposure time for each of the 1200 and 2100 1 



Figure 1. Low-resolution SOAR spectrum from 2013 Jan 15 (0 = 
0.12) showing the late-G/early-K spectral type of the secondary 
and the clear Ho emission. 

mm“^ spectra was 10 minutes. Wavelength calibration 
was performed using FeAr arcs taken after each set of 
three exposures. 

All spectra were reduced in the usual manner, with 
optimal extraction and wavelength calibration using the 
FeAr arcs. To correct for flexure, we cross-correlated 
the sky spectrum for each exposure with a master sky 
spectrum over the wavelength range 6275-6375 A, which 
is rich in sky lines (for the 2100 1 mm“^ data, we used 
the 5577.34 A line). 

We derived a barycentric radial velocity of each spec¬ 
trum through cross-correlation with a set of bright stars 
taken with the same setup, excluding the region around 
Ha in all cases. The median formal uncertainty on the 
radial velocities is 5.2 km s“^. 

2 . 5 . 2 . High-Resolution VLT Spectroscopy 

In cross-correlation of the SOAR spectra with bright 
stars of similar spectra type, we found marginal evidence 
for broadening of the lines, as would be expected if the 
star was rapidly rotating due to tidal locking with the 
primary. However, the resolution of the SOAR spectra 
was too low to precisely measure the line-broadening. 

We obtained a high-resolution spectrum of the optical 
counterpart to IFGL JI4I7.7—4407 on 2014 Aug 20 us¬ 
ing UVES on UT2 of the Very Large Telescope (VLT). 
Three exposures were taken, each 708 sec in length. The 
wavelength coverage was ~ 4800-5750 A and 5850-6800 
A with a resolution of i? ~ 43000. Using the same instru¬ 
mental setup, we also obtained a spectrum of the bright 
K2HI star HD 132096. All spectra were reduced using 
the standard UVES pipeline. To estimate the projected 
rotational velocity (Wot sin i) of the optical counterpart, 
we divided the spectrum of HD 132096 into 100 A chunks, 
convolved these with a set of kernels reflecting a range of 
Vrot sin i values (including limb darkening), and cross- 
correlated these with the unbroadened spectra (Strader 
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et al. 2014). We then fit the relation between the FWHM 
of the cross-correlation peak and the input value of Vrot 
sin i. Finally we cross-correlated the spectrum of the op¬ 
tical counterpart to IFGL J1417.7—4407 with that of the 
standard star in each 100 A region, using the dispersion 
among the measurements as an estimate of the random 
uncertainty in the calculation. 

The value of Vrot sin i we obtain is 33.6 ± 0.7 km s“^. 
We note that this uncertainty only reflects the random 
uncertainty in the calculation, and that the systematic 
uncertainties are likely to be larger. 

3. RESULTS 

3.1. Keplerian Orbit Fitting and Mass Ratio 

As in our previous study of the compact binary 
IFGL J0523.5—2529 (Strader et al. 2014), we performed 
a Keplerian fit to the 47 radial velocities using the IDL 
package BDDTTRAN (Wang et al. 2012) after correcting the 
observation epochs to Barycentric Julian Date (BJD) on 
the Barycentric Dynamical Time (TDB) system (East¬ 
man et al. 2010). Initially we fixed the eccentricity 
e = 0 and fit for the period P, semi-amplitude K 2 , 
systemic velocity Vgys and BJD at superior conjunction 
(To, 5 ; this is when the secondary is behind the primary). 
We found an excellent fit (y^ = 42/47 d.o.f.; rms 5.1 
km s“^), with orbital elements P = 5.37385 ± 0.00035 d; 
K 2 = 115.7 ± 1.1 km s-i; v,ys = -15.3 ± 0.9 km s'^; 
To ,5 = 2457067.605±0.024 d (given as the BJD preceding 
the final spectroscopic dataset). A fit with the eccentric¬ 
ity free yields e = 0 . 01 ± 0 . 01 ,sowe find no significant ev¬ 
idence for a non-zero eccentricity. The spectroscopic and 
photometric periods are consistent to within about 15 s; 
for the remainder of the paper, we adopt the spectro¬ 
scopic period. The phased radial velocity curve is shown 
in Figure 2. 

Using the standard formula for close binaries (Gasares 
2001 ), the mass ratio q = M 2 /Mi is directly determined 
by the rotational velocity and semi-amplitude of the sec¬ 
ondary: Urot sini = 0.462 (1 -I- Using the 

values above, this gives q — 0.179 ± 0.010. 

3.2. Properties of the Secondary 

Given the evidence for an accretion disk and ellipsoidal 
variations, if we assume the secondary is just filling its 
Roche lobe, then the standard formula (Eggleton 1983) 
immediately yields the mean density given the known 
mass ratio and period. This density is p = 0.0068 g 
cm“^. Thus the secondary is clearly a giant. 

The visual impression of the spectrum is that of a late- 
G or early-K spectral type. To quantify this, we used 
the package MKCLASS (Gray & Gorbally 2014), which can 
perform MK classification on flux-calibrated spectra. We 
used two low-resolution spectra taken in 2014 August at 
phases of ^ ^ 0.19 and 0.61. We obtained identical best- 
fit spectral types of G9 for both spectra, and a metallicity 
of [Fe/H] = -0.65. Both of these estimates should be 
taken as uncertain given the low resolution of the spectra. 
Unfortunately, the high-resolution UVES spectrum is not 
suitable for detailed analysis given its modest S/N and 
the rotational broadening. 

We can also study the “night side” properties of the 
secondary with photometry. At this </ = 0, we estimate 
B ~ 17.18 and V ^ 16.22 from the PROMPT photome- 



Figure 2. Phased radial velocities for the optical counterpart to 
IFGL J1417.7—4407, with the listed Keplerian fit overplotted. 

try (see Figure 3; these values are approximate owing 
to the incomplete phase coverage of the light curve). 
Since the system is out of the plane, we assume that the 
full foreground extinction applies. The value given by 
Schlafly & Finkbeiner (2011) is E{B — V) = 0.10. Thus 
we find To = 16.76 and Uq = 15.91, so (T - U)o - 0.85. 
The day side color (</ = 0.5) is {B — V)o 0.90, which 
is slightly redder as expected due to gravity darkening. 



Figure 3. PROMPT photometry in BVR (blue, green, and red 
points) plotted against ellipsoidal models assuming i = 58° and 
Tff = 5000 K. The photometry as plotted is not corrected for 
foreground reddening. 

Using the color-temperature relation of Sekiguchi & 
Fukugita (2000) and a metallicity of [Fe/H] = -0.65, the 
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night side {B — V)o color suggests Teff ^ 5000 K (it would 
be about 200 K lower for solar metallicity). This is con¬ 
sistent with the Tes inferred from a Marigo et al. (2008) 
10 Gyr isochrone of the same metallicity for this {B—V)o. 
Given the unusual evolutionary state of the star with a 
significant amount of mass loss (see below), this value 
should only be taken as a rough estimate. 

3.3. Evidence for Accretion 

As discussed in §2.5.1, IFGL J1417.7—4407 shows per¬ 
sistent double-peaked Ha emission that is evidence for an 
accretion disk in the system. The Ha emission varies in 
its equivalent width and line profile, displaying both sec¬ 
ular and orbital changes. The top panel of Figure 4 shows 
Ha emission line profiles corrected to the rest frame of 
the secondary, averaged over each observing epoch, and 
both normalized and scaled to exhibit variations in the 
profile shape. There are 13 distinct epochs plotted in 
Figure 4 (as discussed in §2.5.1, a subset of the observa¬ 
tions did not include Ha). 

A minimum is consistently seen in the line profile at 
the systemic velocity of the secondary (adjusted here to 
0 km s“^ at all epochs; dashed line in Figure 4). In many 
observations, the emission line is clearly double-peaked— 
a classic signature of an accretion disk. The emission 
appears to be preferentially blueshifted in the first half 
of the orbital period (</> = 0.0-0.5) and redshifted in the 
latter half, although we see a few strong exceptions to 
this trend (e.g., the spectrum from 2013 Jan 15.4). 

The spectrum from 2014 Jun 15.2 {cf = 0.12) shows 
nearly symmetric double-peaked emission, which we take 
as a reference. At this epoch the two peaks are separated 
by 234 km s“^ (marked with dotted lines in Figure 4; 
the blue peak is centered at —112 ± 5 km s“^ and the 
red at 122 ± 5 km s“^). The location of these peaks 
is consistent at many other epochs, although additional 
components with more extreme velocities are sometimes 
superimposed or dominant. 

The observed peak separation of 234 km s“^ is surpris¬ 
ing if we take the simple interpretation of an accretion 
disk such that half the separation (117 km s“^) is the ra¬ 
dial velocity of the outer disk (Smak 1981). Within the 
uncertainties this value is identical to the orbital semi¬ 
amplitude (A '2 = 115 km s“^), implying the disk entirely 
fills the space between the stars. This puzzle is reminis¬ 
cent of that in the black hole binary V404 Gyg, which 
also has a red giant companion. In this system the peak 
separation of the Ha emission components and K 2 are 
similar to those observed in IFGL JI4I7.7—4407. De¬ 
spite intensive study, in V404 Gyg the emission profile 
has not been adequately explained, except to conclude 
that the Ha emission is unlikely to be explained as a 
simple Keplerian disk (Gasares et al. 1993). 

We show variations in the Ha equivalent width as a 
function of orbital phase in the bottom panel of Figure 4. 
Variability is observed on timescales of hours, especially 
on 2015 Feb 18, when the strength and velocity width of 
the emission both ^ double during the night. 

There is evidence of orbital variation, with lower equiv¬ 
alent width around (p ~ 0.5 (compared to ^ 0), but this 
potential signal is contaminated by the observed secular 
changes in the Ha emission. Similar hints of Ha orbital 
modulation are seen in V404 Gyg (Gasares et al. 1993). 


We note that many spectra were not obtained in pho¬ 
tometric conditions and are therefore not flux-calibrated 
in an absolute sense; if the optical continuum varies with 
the Ha emission as might be expected in an accretion 
disk, then equivalent width variations place a lower limit 
on the fluctuation of Ha relative to the companion star 
flux (e.g., Hynes et al. 2002). 



-1500 -1000 -500 0 500 1000 1500 

Ha Radial Velocity (km s ) 



X 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 


Orbital Phase {</>) 

Figure 4. Top: Ho profile as a function of orbital phase (though 
covering many orbits). The UVES spectrum (2014 Aug 20) has 
been smoothed to match the resolution of the SOAR spectra. Bot¬ 
tom: Equivalent width of the Ho emission as a function of orbital 
phase. 

The strongest Ha emission is observed on 2014 May 2 
{(j) = 0.91 and 0.95), accompanied by a redshifted profile. 
These spectra are strongly reminiscent of the redshifted 
Ha “flares” observed in V404 Gyg, which may signify a 
outflow from the accretion disk or central source which 
absorbs the blue wing of the line. In that system these 
variations are hypothesized to arise from magnetic re¬ 
connection events or variable irradiation of the accretion 
disk (Hynes et al. 2002). 

Overall, we interpret the Ha emission as evidence of 
an accretion disk, with rich phenomenology that de¬ 
serves further study. In their “high” low-mass X-ray 
binary states, all three transitional millisecond pulsar 
systems show Ha emission in accretion disks (Wang et 
al. 2009; Pallanca et al. 2013; de Martino et al. 2014). In 
the globular cluster NGG 6397, the companion of PSR 
J1740-5340A shows Ha emission associated with an out¬ 
flow rather than a disk (Sabbi et al. 2003). This star 
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is a stripped giant (Bogdanov et al. 2010; Mucciarelli 
et al. 2013) and so presents an interesting parallel to 
IFGL J1417.7—4407 (see §4) with an evolved compan¬ 
ion but in a rotational-powered rather than accretion- 
powered state. 

3 . 4 . Ellipsoidal Variations and Inclination 

The CSS photometry shows two minima and maxima 
per period, suggesting the presence of ellipsoidal varia¬ 
tions due to tidal deformation of the secondary. Model¬ 
ing the ellipsoidal variations can constrain the inclination 
of the binary, with the largest effect at edge-on inclina¬ 
tions {i = 90°) and no ellipsoidal variations at face-on 
inclinations {i = 0 °). 

However, the very broad CSS filter and moderately 
large uncertainties in the photometry are not ideal for 
modeling the light curve. Therefore, we use the less well- 
sampled but more precise PROMPT photometry in BVR 
for modeling the light curve. 

We make several assumptions in the light curve mod¬ 
eling. First, given the evidence for an accretion disk, 
we assume the secondary fills its Roche lobe. We also 
assume that the accretion disk does not contribute sig¬ 
nificantly to the optical luminosity of the system (this 
is justified by the results below). Due to the modest X- 
ray luminosity of the system (~ 10 ^^ erg s“^; see § 2 . 2 ) 
and long period, irradiation of the secondary is expected 
to be negligible: the X-ray flux at the stellar surface is 
< 10® erg cm“^ s“^, and the ratio of optical to X-ray 
flux at the surface of the star is > 100 ; indeed, we see 
no evidence for irradiation in the light curve. Finally, 
we fix the values of the period and mass ratio to those 
determined via spectroscopy. 

We model the light curve using XRBinary0 (Bayless 
et al. 2010). Civen the assumptions above, the only free 
parameters in our fits are the inclination. Teg, and overall 
normalization of the model. 

Civen the results of §3.2, we begin by fitting models 
with Tes — 5000 K. Photometry from each band was 
modeled independently. The best-fit inclinations using 
the B, V, and R photometry are i = 59 ± 3°, 58 ± 2°, 
and 58 ± 2°. These are clearly in excellent agreement. 
The listed uncertainties in the inclinations are statistical 
only. 

The ellipsoidal models for BVR are plotted against 
the data in Figure 3. We note that the models are not 
perfect representations of the data; for example, at (j) 
between ^ 0 and 0 . 2 , the data appear to slightly lag 
the model, though there is no corresponding feature at 
any other phase. The phase coverage of the data is also 
incomplete, especially in B. The light curve modeling 
should be revisited in the future with improved photo¬ 
metric data. 

To show the effect of changing the temperature, fits 
were also performed with the Teg increased or decreased 
by 250 K. These changes had a minor systematic effect 
on the best-fit inclinations, with i increasing by 1-1.5° 
at Teg = 5250 K and decreasing by the same amount 
for Teg = 4750 K. These give a sense of the systematic 
uncertainties due to Teg on modeling the ellipsoidal vari¬ 
ations. 

® http://www.as.utexas.edu/~elr/XRbinaryV2.4/ 


Another possible systematic error in the photometric 
modeling is the “veiling” of the secondary if there is a 
significant contribution from disk light. Qualitatively, 
veiling leads to an underestimate of the inclination as the 
amplitude of the ellipsoidal modulations is suppressed 
(e.g., Casares et al. 1993). That the same inclination is 
derived for each oi BVR implies either that veiling is 
not important, or that the effective temperature of the 
emitting region of the disk is the same as that of the 
secondary. Future near-IR observations could improve 
constraints on a possible disk contribution to the light 
curve, though veiling may be present even toward redder 
wavelengths in low-mass X-ray binaries (e.g., Kreidberg 
et al. 2012 ). 

For the remainder of the paper we adopt / = 58 ± 2°. 

3 . 5 . Component Masses 

With the measured spectroscopic semi-amplitude and 
period, mass ratio from the rotational broadening, and 
inclination from photometry, we have enough informa¬ 
tion to determine both the primary and secondary masses 
in the system. Our estimate of the primary mass is 
Ml = 1.97 ± O.15M0, consistent with a massive neutron 
star. The secondary is 0.35 ± 0.04Mq. 

We emphasize that the uncertainties on these values 
are solely statistical and do not consider the systematic 
uncertainties described above in the derivation of the in¬ 
clination. However, a substantial increase in the inclina¬ 
tion would be necessary to reduce the inferred primary 
mass to a value consistent with a white dwarf rather than 
a neutron star: i = 72° would imply a primary mass of 
1 . 4 M 0 . Civen the measured mass ratio, the minimum 
primary mass allowed is only 1 . 2 M 0 , so the permitted 
mass range of a putative white dwarf primary is limited. 
In addition, a symbiotic star (accreting white dwarf with 
giant secondary) with a period of 5.4 d would be extraor¬ 
dinarily unusual—most have periods of at least hundreds 
of days (Belczyhski et al. 2000). Future X-ray observa¬ 
tions can provide additional constraints on the inclina¬ 
tion; depending on the geometry of the accretion disk, 
inclinations above i > 70° would produce X-ray dips or 
eclipses. 

3 . 6 . Distance and Kinematics 

The normalization of the best-fit light curve model 
gives the mean bolometric luminosity of the secondary. 
Using the same Marigo et al. (2008) isochrone discussed 
above, the bolometric correction to V is 0.23 mag. Civen 
an observed mean V mag of Vb = 15.80, the implied dis¬ 
tance is 4.4 kpc. The principal uncertainty in this esti¬ 
mate is systematic, and realistically the distance is likely 
uncertain by at least 20 %. 

The Galactic coordinates of IFGL J1417.7—4407 are 
(I ~ 318.9°, b ~ 16.1°): the system is in the direction 
of the inner Galaxy, but above the Plane. The latitude 
corresponds to a height of 1.2 kpc above the Plane for 
a distance of 4.4 kpc. 

Using this distance, the proper motion of 
IFGL J1417.7—4407 (§2.4), and the radial velocity 
derived in §3.1, we can derive the 3-D space motion of 
the system, which is {U,V,W) = (141 ± 74, -137 ± 82, 
-2 ± 96) km s“^ (for this estimate, a Icr distance uncer¬ 
tainty of 0.5 kpc was assumed). While the uncertainties 
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are large due to the corresponding substantial uncer¬ 
tainty in the optical proper motion, the 3-D velocity 
of the binary (219^92^ km s is consistent with that 
observed for other millisecond pulsars (e.g., Gonzalez et 
al. 2011). Given this high velocity, the system’s current 
spatial location may not reflect its origin; nonetheless, 
the modest evidence for subsolar metallicity in §3.2 
would also be more consistent with an origin in the thick 
disk than the thin disk. 

The properties of the system are summarized in Table 

1 . 

3.7. A Chance Alignment? 

Here we briefly consider whether it is plausible that 
the low-mass X-ray binary described in this paper is not 
associated with the Fermi-hAT source that prompted its 
discovery. The total Galactic population of low-mass X- 
ray binaries is estimated to be ~ 10"^ (Jonker et al. 2011), 
with most located in the disk or bulge; of the remainder, 
some will be more distant than the 4.4 kpc estimated for 
this source, or too faint to detect in X-rays. Even op¬ 
timistically assuming a population of ~ 1000 low-mass 
X-ray binaries that meet these constraints, the space den¬ 
sity would only be about 1 per 30 deg^, giving a chance 
alignment probability of < 10“"^ within the 95% Fermi- 
LAT error ellipse of IFGL J1417.7—4407. Even this es¬ 
timate is sanguine, as the X-ray binary is located < 25" 
from the 3FGL FermiAjAA? centroid. We conclude that 
a chance alignment between the low-mass X-ray binary 
and IFGL J1417.7—4407 is very unlikely. 

4. DISCUSSION 

Using spectroscopy and photometry we have shown 
that the Fermi 7 -ray source IFGL J1417.7—4407 is likely 
associated with a low-mass X-ray binary, whose primary 
is likely a neutron star of nearly 2Mq. Given the low 
mass of the hydrogen-rich giant secondary (< O. 4 M 0 ), 
it likely followed a standard evolutionary path for a low- 
mass X-ray binary with a neutron star primary, initiating 
Gase B mass transfer after leaving the main sequence. 
This indicates the system is likely to evolve into a stan¬ 
dard millisecond pulsar with a He white dwarf companion 
at the end of the accretion phase. 

Tauris & Savonije (1999) argue that much of the mass 
in the accretion flow in neutron star binaries with giant 
secondaries is lost rather than accreted onto the neutron 
star. In this case, the inferred high mass of the neutron 
star implies that it was born with a mass above the stan¬ 
dard 1 . 4 M 0 . A similar conclusion has been reached for 
PSR J1614—2230, which has an accurate pulsar mass of 
1.97 ± O.O 4 M 0 measured via the Shapiro delay (Demor- 
est et al. 2010 ), though in this case the secondary is a 
GO white dwarf that likely evolved from an intermediate- 
mass X-ray binary (Tauris et al. 2011). Heavy neutron 
star masses have been inferred for some millisecond pul¬ 
sars discovered via Fermi 7 -ray emission (Romani et 
al. 2012; Schroeder & Halpern 2014). 

Very few low-mass X-ray binaries that are actively ac¬ 
creting have been detected as Fermz-LAT sources (Acero 
et al. 2015), with the exception of the two of the three 
transitional millisecond pulsars (de Martino et al. 2010; 
Stappers et al. 2014). In addition, Bogdanov & Halpern 
(2015) have shown that the Fermi source 3FGL J1544.6- 
1125 is associated with an X-ray binary that exhibits 


optical and X-ray properties similar to the other tran¬ 
sitional millisecond pulsars. The consistent presence of 
7 -ray emission among accreting transitional millisecond 
pulsars motivates us to suggest that IFGL J1417.7—4407 
may also be a transitional millisecond pulsar. There is 
no evidence for a state change in IFGL J1417.7—4407 
over the last ~ 10 years: the GSS photometry shows a 
constant mean magnitude since 2005 (§2.4.1), while PSR 
J1023-I-0038 became brighter by ~ 1 mag during its re¬ 
cent change to the disk state (Halpern et al. 2013). 

The origin of the 7 -ray emission in IFGL J1417.7—4407 
is uncertain. For most millisecond pulsars the 7 -ray emis¬ 
sion is thought to come from the magnetosphere of the 
pulsar and be entirely unrelated to the secondary. The 
lack of radio emission and evidence for an accretion disk 
in this system suggests that this is unlikely to be the cor¬ 
rect explanation in IFGL J1417.7—4407. Indeed, in the 
transitional system PSR J1023-I-0038 the 7 -ray emission 
increased by a factor of at least 5 during the recent state 
change when the pulsar became undetectable (Stappers 
et al. 2014). They hypothesized that the pulsar was en¬ 
shrouded but still active, with the increase in 7 -ray flux 
coming from a shock between the pulsar wind and sur¬ 
rounding material, related to the re-appearance of an 
accretion disk. However, Archibald et al. (2015) show 
that J1023-I-0038 exhibits coherent X-ray pulsations in 
the high state that indicate accreted material is reach¬ 
ing the neutron star. This material ought to quench the 
pulsar, so it is not trivial to associate the 7 -ray emission 
with a pulsar wind. An alternative model is that the 
7 -rays are due to the interaction between a relativistic 
jet and the surrounding material, or to synchrotron self- 
Gompton radiation produced at the boundary of the in¬ 
ner accretion disk (Papitto et al. 2014; Deller et al. 2015). 

The X-ray luminosity of IFGL J1417.7—4407 (~ 1.4 x 
1033 gj.g g-i^ jg consistent with the average value ob¬ 
served in the disk state of the known transitional millisec¬ 
ond pulsars (Linares 2014), and could well suggest that 
mode switching is also occurring in IFGL J1417.7—4407, 
but is not resolved in our short X-ray observations due to 
the greater distance and hence lower flux of this source. 
The X-ray photon index (P ~ 1.3) is also consistent with 
the transitional objects in this state. The ratio of 7 -ray 
to X-ray luminosity is ~ 20, which is substantially larger 
than that observed in the disk state for J1023-I-0038, per¬ 
haps related to the giant companion. 

The transitional millisecond pulsars, including 
J1023-I-0038 (Deller et al. 2015), have been detected 
as flat-spectrum cm radio continuum sources in the 
low-mass X-ray binary state. The observed 10 GHz 
radio continuum flux density of J1023-I-0038 would 
correspond to ~ 5-10 /iJy at the inferred distance of 
IFGL J1417.7—4407, which is marginally detectable 
with southern radio interferometers. 

Glearly, additional work is needed to characterize 
IFGL J1417.7—4407. Gontinuing observations with 
Fermz-LAT could allow a measurement of 7 -ray variabil¬ 
ity of the source, which would help determine the ori¬ 
gin of the 7 -ray emission. Radio observations at superior 
conjunction should be used to search for radio pulsations 
or provide strong constraints on their presence, although 
radio pulsar emission is not expected while the source 
has an accretion disk. Improved photometric data and 
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Table 1 

Summary of Properties 


Opt. R.A. (J2000 h:m:s) 
Opt. Dec. (J2000 ° 
Period (d) 

K 2 (km s“^) 

To. 5 (d) 

Vrot sin i (km s 
M 2 1 Ml 

in 

a (Rq) 

Vo, 4> = 0 (mag) 

{B - V)o, <j) = 0 (mag) 
Ml (Mq) 

M2 (Mq) 

Vaya (km S“l) 
fiaCosS (mas yr“^ ) 
fis (mas yr-l ) 

(U,V,W) (km s-i) 
Distance (kpc) 


14:17:30.60 
-44:02:57.4 
5.37385 ± 0.00035 
115.7 ± 1.1 
2457067.605 ± 0.024 
33.6 ±0.7 
0.179 ±0.010 
58 ±2 
17.1 ±0.4 
15.91 
0.85 
1.97 ±0.15 
0.35 ±0.04 
-15.3 ±0.9 
-8.8 ±5.0 
-3.1 ±4.7 

(141 ± 74, -137 ± 82, -2 ± 96) 
4.4 


modeling could aid in constraining the systematic uncer¬ 
tainties on the mass estimate of the neutron star. Ob¬ 
servations of X-ray variability on short timescales, an 
improved X-ray spectrum, and the detection of radio 
continuum emission offer the best current hope for de¬ 
termining whether IFGL J1417.7—4407 is a transitional 
millisecond pulsar. 
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